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SUMMARY 


The  variational  principle,  differential  equations,  and  boundary 
conditions  governing  the  cross-sectional  distortions  due  to  inertia 
loading  of  a  two-dimensional  model  of  a  thin  monocoque  wing  are  shown. 
A  theoretical  analysis  of  this  simplified  model  is  made  in  order  to 
determine  the  nature  of  the  coupling  between  the  cross-sectional  modes 
and  the  spanwlse  deformation  modes.  General  solutions  are  obtained  in 
finite-difference  form  for  arbitrary  cross  sections  and  an  exact  solu¬ 
tion  is  presented  for  a  parabolic-arc  cross  section  of  constant  cover 
thickness.  The  application  of  these  results  in  evaluating  the  coupled 
frequencies  of  the  act\ial  structure  is  discussed.  Frequencies  evalu¬ 
ated  for  a  parabolic-arc  monocoque  beam  show  good  agreement  with 
experimental  values. 


INTRODUCTION 


Chordwise  modes  -  that  is,  modes  whose  chordwise  variation  in 
deflection  is  such  that  more  than  one  node  line  exists  at  every  cross 
section  along  the  span  -  have  been  observed  in  several  dynamic  inves¬ 
tigations  of  pure  monocoque  and  multiweb  wing  models.  For  example,  in 
an  extensive  series  of  tests  of  aero dynamically  heated  multi web  wing 
structures  the  principal  mode  of  failure  was  established  to  be  chord- 
wise  flutter.  (See  ref.  1.)  Also,  in  an  investigation  of  the  vibra¬ 
tional  characteristics  of  circular-arc  monocoque  beams,  discussed  in 
reference  2,  several  "extra"  resonant  conditions  were  observed  experi¬ 
mentally  which,  although  referred  to  as  panel  modes,  are  more  accu¬ 
rately  described  as  chordwise  or  cross-sectional  modes.  In  addition, 
the  results  of  reference  2  indicate  that  the  accuracy  of  the  predicted 
overall  bending  and  torsion  frequencies,  even  after  the  inclusion  of 
inp)ortant  secondaiy  effects,  such  as  transverse  shear  and  shear  lag 
in  bending  and  restraint  of  warping  in  torsion,  are  considerably 
impaired  by  the  presence  of  these  chordwise  modes.  For  instance,  the 


calculated  frequency  of  even  the  first  syiimetrical  torsion  mode  (mode  2 
in  table  II  of  ref.  2)  differed  from  the  one  obtained  exq)erimentally  by 
72  percent  while  the  third  symmetrical  mode  differed  by  more  than 
250  percent. 

The  influence  of  cross-sectional  distortions  on  the  overall  beam 
modes  was  investigated  in  references  3  and  4.  These  investigations^ 
however^  were  not  concerned  primarily  with  the  prediction  of  cross- 
sectional  modes  but  rather  with  the  influence  of  cross-sectional  flexi¬ 
bilities  on  the  overall  beam  modes.  Furthermore^  neither  analysis  is 
applicable  to  monocoque  or  multlweb  wings  with  thin  airfoil-type  cross 
sections  such  as  those  used  in  references  1  and  2. 

The  purpose  of  the  present  paper  is  to  investigate  the  chordwise 
modes  and  frequencies  of  beams  with  thin  airfoil-type  cross  sections 
and  to  determine  the  influence  of  cross-sectional  flexibilities  on  the 
overall  beam  frequencies.  A  sinplified  model  which  embodies  the  prin¬ 
cipal  mechanisms  of  cross-sectional  and  spanwlse  flexibility  is  analyzed. 
The  differential  equations  and  boundary  conditions  for  the  simplified 
model  appropriate  for  wings  of  arbitrary  airfoil-type  cross  section  are 
derived.  An  exact  solution  is  foimd  for  the  model  for  the  particular 
case  of  a  monocoque  wing  of  parabolic  cross  section.  In  addition,  a 
general  solution  is  presented  in  finite-difference  form  by  applying  the 
calculus  of  variations  directly  to  the  energy  expressions.  A  discussion 
is  made  of  the  use  of  these  results  to  predict  the  frequencies  of  both 
the  coupled  cross-sectional  modes  and  the  overall  beam  modes  of  a  mono¬ 
coque  wing  presented  in  reference  2. 


SYMBOI£ 


c 

E 

h 

tio 

I 

ki,k2,k 

M 


chord  length  of  beam 
modulus  of  elasticity 

distance  from  x-axis  to  median  line  of  skin 
maximimi  value  of  h 

moment  of  inertia  about  median  line  of  skin,  t5/l2 
spring  constants 

moment  (see  fig.  2) 


N 


niunber  of  discrete  stations  in  finite-difference  analysis 


radius  of  curvature  of  midplane  of  beam 
tangential  coordinate 
maximum  kinetic  energy- 
skin  thickness 
maximum  strain  energy- 

displacement  of  midplane  in  x-direction 
displacement  in  s-direction,  h'w  +  u 
shear  (see  fig.  2) 

deflection  in  z-dlrection  (positive  upward) 
deflection  normal  to  s-direction 
chordwlse  coordinate 
spanwise  coordinate 

coordinate  perpen'ilcular  to  x  and  y  coordinates 
frequency  coefficient 

frequency  coefficient  for  uncoupled  chordwlse  mode 
strain  in  s-directlon 
strain  in  x-dlrectlon 

nondlmensional  distance  bet-ween  stations  on  the  half-chord 
normal  midplane  force  in  the  beam 

1 

mass  density  of  cover 
natural  frequency 

uncoupled  cross-sectional  frequency 


“c,l 

^exper 


fundamental  or  first  uncoupled  cross-sectional  frequency 

experimental  natuiral  frequency 

rigid-tody  bending  frequency,  \ 

y  cpt 


^r, T  rigid-body  torsion  frequency. 

Primes  on  symbols  denote  differentiation  with  respect  to  x;  primes 
on  matrices  denote  the  transpose  of  a  matrix. 


THEORETICAL  ANALYSIS 


Idealization 

A  rigorous  approach  to  the  problem  of  the  vibrational  characteris¬ 
tics  of  a  monocoque  wing  which  includes  both  spanwise  and  cross-sectional 
flexibilities,  would,  of  necessity,  involve  some  type  of  shell  analysis. 

A  more  appealing  euLthough  less  rigorous  approach,  however,  is  to  investi¬ 
gate  first  the  uncoupled  chordwise  and  spanwise  characteristics  of  the 
structure  independently  and  then  determine  the  influence  on  the  natural 
frequencies  of  the  coupling  of  these  motions.  This  can  be  done  in  a 
manner  similar  to  that  used  in  references  5  and  4-  by  mounting  a  typical 
cross  section  of  the  wing  on  elastic  supports  in  such  a  way  as  to  allow 
for  the  necessary  overall  degrees  of  freedom. 

The  only  structiires  considered  in  this  paper  have  uniform  properties 
in  the  spanwise  direction  with  a  symmetrical  curved  airfoil  cross  section, 
as  shown  in  figxire  l(a)  .  The  wing  or  beam  is  of  either  pure  monocoque  or 
raultiweb  construction  with  no,  or  relatively  few,  bulkheads.  An  Idealized 
structure  that  can  be  easily  analyzed  but  still  possesses  essentially  the 
same  cross-sectional  and  spanwise  flexibilities  as  the  original  is  shown 
in  figiire  l(b) . 

The  idealized  structure  consists  of  two  curved  beams  attached  at 
the  ends  to  elastic  supports  in  such  a  manner  as  to  maintain  the  angle 
between  the  beams  at  each  end.  The  idealization  is  the  same  for  both 
the  monocoque  and  multiweb  prototypes.  Consequently,  the  spanwise  webs 
of  the  multiweb  structiire  are  assumed  to  contribute  nothing  to  the  chord- 
wise  stiffness  of  the  wing.  The  webs,  however,  do  tend  to  prevent  any 
relative  motion  between  the  top  and  bottom  covers.  Hence  the  additional 
assumption  is  made  that  deflections  of  the  top  and  bottom  beams  in  the 
idealized  structure  are  the  same.  The  depth  of  the  individual  beams  in 


_f_  J_, 
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the  idealized  structure  is  taken  to  he  the  same  as  the  thickness  of  the 
covers  in  the  prototype.  The  stiffness  of  the  elastic  supports  or  springs 
is  chosen  to  duplicate  the  bending  or  torsional  stiffness  of  the  proto¬ 
type,  depending  on  whether  the  motion  is  symmetrical  or  antisymmetrical 
about  the  midchord. 


L 


Method  of  Analysis 

The  structure  to  be  analyzed  is  shown  in  figure  2  along  with  the 
coordinate  system  and  sign  convention  to  be  used.  As  the  structure  is 
symmetrical  about  the  x-axis  and  the  distance  between  the  top  and  bottom 
segments  is  maintained,  only  half  of  the  structure  need  be  considered. 
Each  segment  is  assumed  to  carry  midplane  forces  in  addition  to  trans¬ 
verse  shear  and  bending  moments  (see  fig.  2),  and  will  be  considered  to 
act  in  a  manner  consistent  with  elementary  beam  theory  (i.e.,  plane  sec¬ 
tions  remain  plane). 

The  strain  in  the  midplane  of  the  beam  is  given  by 


®s  " 


du 

ds 


(1) 


where 

w  deflection  normal  to  s-direction 

r  radius  of  curvature  of  midplane 

u  displacement  in  s-direction 

s  tangential  coordinate 

Since  the  cross  section  of  the  prototype  being  analyzed  is  that  of 
a  thin  airfoil  section  and  deformations  are  considered  small,  the  tan¬ 
gential  coordinate  s  can  be  approximated  by  the  chordwise  coordinate  x, 
and  the  curvature  by  -h" .  Furthermore,  the  deflection  w  can  be  con¬ 
sidered  to  be  the  same  as  the  deflection  in  the  z-direction,  and  the  dis¬ 
placement  in  the  s- direction  can  be  written  as 

u  =  h'w  +  u  (2) 

where  w  is  the  deflection  in  the  z- direction  and  u  is  the  displace¬ 
ment  in  the  x-direccion.  Hence  the  expression  for  strain  in  equation  (l) 
becomes 


e„  =  h’w'  +  u' 


(3) 
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If  this  expression  is  used  for  the  strain,  and  w  and  u  are  considered 
as  amplitudes  of  displacement,  the  maximijm  strain  energy  for  one  of  the 
covers  of  the  idealized  structure  vibrating  in  a  natural  mode  becomes 


2  0 

f  El(w")^dx  +  f  A(h'w'  +  u' 

'0  *^0 

)dx  +  i  kiw2(c)  +  i  k2W^(0)  (4) 

where 

El 

bending  stiffness  of  beam 

L 

1 

A 

normal  midplane  force 

4 

4 

ki,k2 

spring  constants 

4 

The  first  term  in  equation  (4)  represents  the  contribution  to  the  strain 
energy  due  to  bending  of  the  idealized  structure,  the  second  represents 
the  strain  energy  due  to  the  midplane  force,  and  the  third  and  fourth 
terms  are  the  potential  energy  of  the  spring  mounts. 

In  the  present  analysis  the  force  A  will  be  considered  constant 
and  equal  to  the  average  midplane  force;  hence,  with  the  use  of  equa¬ 
tion  (3), 


A  = 


+  u' )dx 


(5) 


Furthermore,  if  inextensional  theory  is  assumed  u'  vanishes  in 
equation  (4),  and  equation  (5)  becomes 


A 


1 

c 


Eth’w'dx 


(6) 


The  majcimum  kinetic  energy  of  the  idealized  structure  vibrating  in 
a  natural  mode  is 


(7) 


0 


where  cn 
tlon  and 


is  the  natural  circular  frequency  of  the  mode  under  considera- 
p  is  the  mass  density  of  the  cover. 
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A  natural  mode  of  vibration  must  satisfy  the  variational  equation 

6(U  -  T)  =  0  (8) 

where  the  variation  is  taken  independently  with  respect  to  w.  Applica¬ 
tion  of  equation  (8)  and  the  usual  methods  of  the  calculus  of  variations 
results  in  the  following  differential  equation  and  boundary  conditions: 

(9) 
(10) 

(11) 

(12) 

(15) 

where  A  is  given  by  equation  (6) . 

If  the  cover  thicknesses  of  the  original  wing  are  constant  and  if 
the  half-depth  h  of  the  wing  is  approximately  parabolic  (as,  for 
example,  the  circular-arc  airfoils  used  in  refs.  1  and  2)  an  exact  solu¬ 
tion  to  the  differential  equation  is  possible  and  is  discussed  in  the 
next  section.  However,  for  an  airfoil  of  arbitrary  shape  and  variable 
cover  thickness  it  is  unlikely  that  an  exact  solution  can  be  found. 

For  these  cases  the  use  of  some  approximate  variational  procedure  is 
advisable.  One  such  approach  is  presented  in  the  appendix.  Here,  by  a 
procedure  similar  to  that  discussed  in  reference  the  energy  is 
expressed  in  matrix  form  in  terms  of  the  deflections  at  discrete  points 
by  means  of  finite- difference  relationships.  The  variational  procedure 
results  in  a  frequency  equation  in  matrix  form  from  which  the  mode 
shapes  and  frequencies  of  the  structure  can  be  calculated. 


(EIw")"  -  Ah"  -  ptcn^w  =  0 


(EIw")5w' 


=  0 


x=c 


(EIw")’  -  Ah'  -  w 


6w 


=  0 


x=c 


(EIw")5w'  =  0 

x=0 


(EIw")'  -  Ah'  w 


6w 


x=0 


=  0 


Exact  Solution  of  Idealization  of  the  Monocoque  Be.am 
of  Parabolic-Arc  Cross  Section 

For  an  airfoil  of  doubly  symmetrical  cross  section  and  constant 
cover  thickness  whose  half-depth  h  can  be  represented  by  a  parabolic 
equation,  an  exact  solution  to  the  differential  equation  for  the  cross 
section  mo\mted  on  springs  is  possible.  In  order  to  take  advantage  of 


8 


( 


the  double  symmetry  the  origin  of  the  coordinate  system  is  now  relocated 
at  the  midchord  of  the  cross  section  and  the  nondimensional  coordinate 

Py 

I  =  -^ —  1  is  used.  Thus,  the  half-depth  can  be  represented  as 


h  =  ho(l  -  (Ik) 

Furthermore,  the  spring  constants  are  identical,  so  that  =  k2  =  k. 
Equation  (9)  can  now  be  written  in  nondimensional  form  as 


where 


d^w  Ac^  d^h  k  n 

— -  -  1 - a  w  =  0 

d|4  kEl 


k  2 
Ij.  c  ptoi 

16EI 


(15) 


(16) 


A  - 


4hotE 


w(-l)  -  w(l)  + 


w(|)d| 


-1 


(17) 


Because  of  the  constant  cover  thicknesses  of  the  original  wing,  the 
more  stringent  conditions  of  inextensional  theory  are  not  necessary  in 
determining  A.  It  need  only  be  assumed  that  there  is  no  overall  short¬ 
ening  or  stretching  of  the  cross  section;  then  the  integral  of  u'  will 
vanish  in  equations  (4)  and  (5)  instead  of  u'  being  zero  everywhere. 


Because,  in  the  general  solution  to  the  differential  equation,  the 
symmetrical  and  anti symmetrical  parts  are  easily  recognizable,  they  are 
presented  separately. 


Symmetrical  solution.-  For  the  symmetrical  solution  the  boundary 
conditions,  as  given  by  equations  (lO)  and  (ll),  can  be  written  in  non- 
dimensional  form  as 


and 


d^ 

d|2 


1=1 


=  0 


(18) 


=  0 


L 

1 

4 

4 

4 


1=1 


(19) 
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The  synunetrical  solution  to  the  homogeneous  part  of  equation  (l5)  is 

w  =  cos  04  +  Bp  cosh  a4  (20) 

emd  the  particular  solution  is  simply 


w  = 


AhoC 

2EIa^ 


(21) 


Thus 


w  =  Ap  cos  a5  +  Bp  cosh  + 


TlioC^ 


2EIa^ 


(22) 


Substitution  of  equation  (22)  into  the  boundary  conditions  for  the  sym¬ 
metrical  solution  (eqs.  (l8)  and  (l9))  yields,  respectively, 


.  cosh  a  _ 

Ap - Bp 

cos  a 


(25) 


Bl 

cosh  a  = 
A 


_  i  ^  _ 

^  \  l+Et^a^, 


QaPf 


El  ^::~(tan  a  +  tanh  a)  -  k 
,5 


i2k) 


Note  that  three  unknowns  appear  in  equations  (25)  and  (24):  Ai, 
Bp,  and  A.  However,  by  substituting  equation  (22)  for  w  into  equa¬ 
tion  (17)  for  A  and  performing  the  necessary  integration,  a  third 
expression  is  obtained: 


8h  tE  / 

A  =  — ^ —  Bp  cosh  a^-2  + 


tan  g  +  tanh 
a 


-) 


(25) 


If  Bp  in  equation  (25)  is  replaced  by  its  identity  in  equa- 

A 

tlon  (24),  the  characteristic  or  frequency  equation  is  found,  and  can 
be  written  as 


1  = 


(26) 
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Substitution  of  equations  (25)  and  (25)  into  equation  (22)  yields  the 
mode  shape 


w 


cos  a|  cosh  04 
cos  a  cosh  a 


tan  a  +  tanh  a 

5k  /c\5 

a 

2Eat\t/ 

1  5k  /c 

\5 

4-Ea^  v; 

) 

(27) 


Antisymmetrical  solution.-  In  the  antisymmetrical  solution  the  mid¬ 
plane  force  A  is  zero,  as  can  be  seen  from  equation  (6) .  In  equation  (6) 
the  integrand  of  the  integral  expression  for  A  is  an  antisymmetrical 
function  about  c/2;  consequently  its  integral  from  0  to  c  is  zero. 

The  differential  equation  for  the  cross  section,  given  by  equation  (l5)> 
and  the  boundary  conditions  given  by  equations  (lO)  and  (ll)  therefore 
reduce  to 


and 


k 

a  w 


0 


(28) 

(29) 


d^v 

d|5 


(50) 


Equations  (28)  to  (jO)  a^re  the  differential  equation  and  bovindary 
conditions  for  a  beam  supported  by  springs  at  |  =  -1  and  |  =1.  The 
solution  is  presented  in  most  vibration  text  books,  but  for  the  sake  of 
completeness  it  will  be  repeated  here.  For  the  antisymmetrical  solution 


w  =  A2  sin  04  +  B2  sinh  04  (5^) 

Substitution  of  equation  (51)  into  the  boundary  conditions  yields 


A2 


sinh  a  ^ 


(52) 


cos  a  + 


B2 


cosh  a)  =  ^{^2  a,  +  B2 


sinh  a) 


(55) 
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Thus,  the  characteristic  or  frequei^cy  equation  becomes 

coth  a  -  cot  a  _  /c\5 


and  the  mode  shape  is 


sin  cxf  sinh  of 

- -  +  - - 

sin  a  sinh  a 


(5^) 


(55) 


NUMERICAL  RESULTS  FOR  IDEALIZATION  OF  MONOCOQUE  WING 
WITH  SYMMETRICAL  PARABOLIC  CROSS  SECTION 


Symmetrical  Modes 


The  exact  frequency  and  mode- shape  equations  of  the  idealized  struc¬ 
ture  were  obtained  in  the  previous  section.  In  equation  (26)  the  fre- 


a 


quency  parameter 

\ 

given  as  a  fmction  of 


where 


(ii  =  a~||~u~~  I  of  the  symmetrical  modes  was 

t/hg  and  the  parameter  .  The  motion 

of  the  idealized  structure  can  be  considered  to  be  the  result  of  the 
coupling  of  two  types  of  modes  of  vibration:  an  infinite  set  of 
uncoupled  cross-sectional  modes  which  are  obtained  by  omitting  the 
springs,  and  a  rigid-body  mode  represented  by  a  rigid  cross  section 
mounted  on  springs.  Consequently,  if  the  vibrations  are  considered  in 

this  light  a  more  useful  parameter  than  wo-uld  be 

where  cUj.  g  is  the  rigid-body  translation  frequency  and  is  defined 

so  that 


,  0)  ) 


Thus  the  frequency  equation  (26)  becomes 


48 


1 


_ 2a 

tan  a  +  tanh  a 


+  1 


1 


_ 2a _ 

tan  a  +  tanh  a 


(56) 


(57) 
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and  the  corresponding  mode  shape  (eq.  (27))  is  given  by- 


cos  a4  cosh 

w  =  - ^  f  - ^ 

cos  a  cosh  a 


a  +  tanh  a 

a  "  CD 

.  -  (^f 

V  (’ll  / 


The  coupled  frequency  coefficient  a  was  calculated  from  equa¬ 


tion  (57)  for  various  values  of 


and  the  results 


are  shown  in  figure  3*  Here  the  frequency  coefficient  a  for  the  first 

/  CD  / 1 

three  modes  is  plotted  as  a  function  of  / - )  for  a  range  of  (^) 


from  0  to  0 . 09 .  As 


becomes  greater  than  1,  all  curves  approach 


horizontal  asymptotes  whose  values  depend  on  ^ —  These  asymptotes  have 

ho 

the  value  of  the  uncoupled  cross-sectional  frequencies  obtained  by 
letting  k  =  0  which  correspond  to  the  case  where  no  coupling  exists 
between  rigid  and  cross-sectional  motion.  Consequently  the  mode  shapes 
associated  with  these  frequencies  are  ones  which  are  predominantly  cross 
sectional  in  character  -  that  is,  they  contain  two  or  more  nodal  points. 

For  k  equal  to  zero,  equation  (37)  reduces  to  the  frequency  equa¬ 
tion  for  the  uncoupled  cross-sectional  mode: 


\hQ  /  ^4\tan  ttg  +  tanh  j 

and  equation  (38)  reduces  to  the  corresponding  uncoupled  cross-sectional 
mode  shai)e: 

cos  a-|  cosh  a„|  2  . 

V  =  ■  ■  ^  ^  I  ■■  .  ■  .  ■  I  -  - .  ( 40 ) 

cos  cosh  Op  oc^  / 


where 


L 

1 

k 

k 

4 
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In  figures  4(a)  and  4(b),  Is  shown  plotted  as  a  function  of 

for  the  first  two  symmetrical  frequencies.  These  plots  serve  to 

determine  the  uncoupled  chordwise  frequency  which  corresponds  to  the 
pure  chordwise  modes  of  an  infinitely  long  wing. 


Values  at 


CJD 


03. 


’r,B 


=  0  in  figure  3  describe  the  condition  where 


the  ends  of  the  cross  section  are  pinned;  that  is,  k  =  oo.  Hence,  mode 


shapes  for 


03 


03. 


'r,B 


«  1  are  also  of  predominantly  chordwise  character. 


For  values  of 


03 


03. 


'r,B 


approximately  equal  to  1  the  curves  approach  a 


vertical  asymptote  =  1  and  correspond  to  the  condition  where 

\“r,  B/ 

the  mode  shapes  are  of  the  rigid-body  type;  that  is,  they  contain  no 


nodal  points.  As  can  be  seen  from  figure  5;  when 


=  0  the  curves 


fall  on  the  asymptotes  and  no  coupling  exists.  This  uncoupling  of  the 
two  motions  could  also  have  been  surmised  by  examining  the  mode  shape 

^2 


as  given  by  equation  (38)-  Note  that  for 


_t_ 


=  0  the  tip  deflec¬ 


tions  of  the  cross  section  in  equation  (58)  are  zero,  and  hence  the 
stiffness  of  the  springs  can  have  no  influence  on  the  overall  motion. 
The  differences  between  the  curves  and  their  corresponding  asymptotes 

for  the  other  values  of  (  — )  represent  the  change  in  frequency  due 

\^/ 

to  coupling  between  the  chordwise  and  rigid  motions.  As  can  be  seen 
from  the  plots,  the  greatest  effect  is  felt  when  U3  < 


One  drawback  with  the  curves  presented  in  figure  5  Is  the  difficulty 
encountered  in  using  these  curves  to  calculate  the  coupled  frequency, 
since  the  iinknown  coupled  frequency  coefficient  a  appears  in  both  the 
abscissa  and  ordinate  of  the  curves.  Thus  the  results  of  figure  3  have 
been  replotted  in  figures  5(a.)  and  5(1')  in  terms  of  slightly  different 
parameters. 

In  figijires  5(a)  and  5(l)  He  ratio  of  the  coupled  frequency  cd  to 
the  rigid-body  frequency  cu^  g  is  plotted  as  a  fxxnction  of  the  ratio 

of  the  fundamental  uncoupled  chordwise  frequency  03  to  the  rigid- 

± 

“c  1 

body  frequency  a3„  .q.  Figure  5(a)  presents  the  data  for  values  of  — - — 
r#  n  a3„  ^ 
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from  0  to  4.8  and 


O). 


’r,B 


from  0  to  3-2.  Figure  5(b)  is  simply  an 


1 

enlargement  of  figure  5(a)  in  the  region  where  — —  goes  from  0  to  2.0. 

“r,B 

The  values  of  in  the  abscissa  are  based  on  the  first  uncoupled 

symmetrical  mode  for  k  =  0.  Hence  the  abscissa  depends  only  on  the 
stiffness  properties  of  the  structure  and  not  on  the  frequency  parameter. 
The  vertical  asymptote  of  figure  3  appears  in  figure  5  as  a  horizontal 


line  at 


0) 


CD 


=  1.0,  while  the  horizontal  asymptotes  of  figure  3  trans- 


r,B 


form  into  straight  lines  which  emanate  from  the  origin.  Note  that  only 
/t 

the  (  — )  =0  asymptotes  appear  in  figure  5>  since  the  asymptotes  for 

\bo/ 

2  2 
all  other  values  of  ( — )  lie  in  such  close  proximity  to  the  fr— )  =  0 

\^o)  \^o/ 

asymptotes  that  they  cannot  be  easily  distinguished. 

Figure  5  shows  that  for  values  of  n  greater  than  co  the 

influence  of  coupling  is  appreciable  only  for  the  first  or  lowest  mode, 


while  for 


Ob. 


'c,l 

'r,B 


<  1  all  modes  investigated  show  substantial  effects. 


The  frequencies  of  the  coupled  modes  are  easily  obtained  from  figure  5; 
since  the  abscissa  depends  only  on  the  physical  properties  which,  of 
course,  are  assumed  to  be  known.  The  difficulty  lies  in  establishing 
the  characteristics  of  the  mode  shapes  -  that  is,  in  determining  whether 
the  resulting  mode  is  predominantly  a  rigid  mode  with  no  nodal  points  or 
predominantly  a  chordwise  mode  with  two  or  more  nodal  points. 


In  order  to  try  to  resolve  this  problem,  the  first  three  modes  for 


^  *1 

-i-')  =  0.018  and  =  0.987  and  0.409  were  calculated  and  are 

^0/  ^r, B 

shown  in  figures  6(a)  and  6(b).  From  figure  6(a)  it  can  be  seen  that 
the  first  mode  is  predominantly  rigid-body  motion  while  the  remaining 
modes  consist  of  chordwise-type  motions.  In  figure  6(b),  however,  the 


second  mode  is  the  rigid-body- type  mode.  Presumably,  as 


CD 


cj. 


%B 


is 


reduced  still  further  the  third  mode  becomes  the  rigid-body  mode  and  so 
on.  In  order  to  illustrate  the  change  in  mode  shape  due  to  coupling  of 
chordwise  and  rigid-body  motions,  a  sketch  of  the  second  mode  for 


“c  1 

A  =  0.018  o-iu  various  values  of  ' 


W) 

is  shown  in  figure  6(c). 


CD 


ranging  from  0.026  to  O.987 


r,B 


4=-  4=-  4=-  H  f 


4^  4=-  H  f 
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From  the  foregoing  discussion  it  can  be  seen  that  for 


CD 


CD. 


S-2-i  >  1  the 


'r,B 


CD, 


rigid  mode  is  alvays  the  lowest  mode.  For 


'c,l 


CD 


<  1  it  is  advisable  to 


r,B 


calculate  the  mode  shape  in  order  to  determine  which  mode  is  actually 
the  rigid-body  mode. 


In  reference  k  the  behavior  of  rectangular  bents  with  flexible 
members,  mounted  on  either  deflectional  or  rotational  springs,  has  been 
analyzed  in  detail. 


It  is  of  interest  to  note  that  the  assumption  of  hQ  =  0  =  ooj 

in  the  present  paper  corresponds  to  the  limiting  case  in  reference  4 
where  the  depth  of  the  web  is  assumed  to  become  zero  while  the  thickness 
becomes  infinite.  For  such  a  configuration  the  frequency  equations  of 
the  two  papers  are  found  to  be  identical.  This  is  true  for  both  the 
symmetrical  and  the  antisymmetrical  modes. 


Antisymmetrical  Modes 


For  the  antisymmetrical  modes,  the  frequency  parameter  a  as  given 

in  equation  (54)  is  a  function  of  only.  This  lack  of  dependence 

of  frequency  on  camber  in  the  antisymmetrical  modes  is  a  consequence  of 
the  vanishing  of  the  midplane  force  A.  Therefore,  the  equations  for 
the  mode  shapes  and  frequencies  are  those  of  the  antisymmetrical  vibra¬ 
tions  of  a  beam  mounted  on  spring  supports.  As  was  true  for  the  symmet¬ 
rical  modes,  the  motion  of  the  idealized  structure  can  be  considered  as 
a  result  of  the  coupling  between  cross-sectional  modes  and  rigid-body 
modes.  In  the  antisymmetrical  case,  however,  the  rigid-body  mode  is  a 


pitching  or  torsional  mode. 
f^T  T^ 

is  I J  ,  in  which  ^ 


Hence,  a  more  useful  parameter  than 
is  now  defined  by 


CO. 


r,T 


5k 

ptc 


(42) 


The  parameter  cjd  „  can  be  looked  upon  as  a  rigid-body  torsion  mode 
1 

represented  by  the  well-known  equation 


l6 


CD. 


I. 


(^5) 


where 


and 


K 


ptc' 


(H) 


(45) 


Note  that  Ip  is  the  approximate  polar  moment  of  inertia  for  a  thin 

cross  section.  The  characteristic  or  frequency  equation  (eq.  (54))  can 
then  he  written  as 


3  coth  a  -  cot  a 
2  a 


CD. 


'r,  T 

O) 


(46) 


The  mode  shape  is  not  an  explicit  function  of  the  spring  constant  and 
remains 


w  = 


sin  a|  slnh  a| 
sin  a  sinh  a 


(47) 


The  coupled  frequency  coefficient  a  was  calculated  from  equa- 

2 

tion  (46)  for  various  values  of  ( - 1  and  the  results  are  shown  in 

\“r,T/ 

figure  7*  In  figure  7  the  first  three  modes  of  a  are  plotted  as  a 

2 

function  of  (  — )  .  The  uncoupled  cross-sectional  modes  are  repre- 

sented  hy  dashed  horizontad  lines,  at  a  =  5*927  and  7*069,  for  k  =  0. 
The  pin-ended  condition  for  k  =  oo  is  represented  hy  the  dashed  lines 
at  a  =  rt,  2it,  and  3^*  The  rigid-hody  torsion  mode  is  the  dashed  vertl- 

^  =1.  The  coupled  modes  are  shown  hy  the  curved  lines. 


cal  line  at 


0). 


'r,  T 


The  corresponding  plot  of 
figure  8,  The  values  of 


CD 


CD 


%T 


as  a  function  of 


c,l 

“V,T 


is  shown  in 


CD. 


'c,l 


in  the  abscissa  are  based  on  the  first 


uncoupled  antlsymmetrical  mode  for  k  =  0.  The  results  are  plotted  in 
a  manner  similar  to  that  for  the  symmetrical  modes  hut  only  one  curve 
is  presented  for  each  mode,  the  results  being  independent  of  the  ^ 


or  camber  effect.  The  dashed  horizontal  line  at 


CD 


av,T 


=1.0  represents 


-P"  -P"  -fr  M  tr* 


F 


IT 


the  rigid-body  torsion  mode,  and  the  dashed  strai^t  lines  inclined  to 
the  abscissa  represent  the  uncoupled  chordwise  modes. 

Contrary  to  the  case  of  the  symmetrical  modes,  the  predominantly 
rigid-body  mode  with  no  nodal  points  is  easily  distinguished,  as  it  is 
always  the  first  or  fiindamental  mode. 


L 

1 

4 

4 
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APPLICATION  OF  METHOD 


Ihe  present  paper  has  two  objectives:  (l)  to  analyze  the  chordwise 
vibrational  characteristics  of  thin  wings  with  airfoil-type  cross  sec¬ 
tions  and  (2)  to  obtain  the  effect  of  the  coupling  between  the  chordwise 
and  spanwlse  motions  of  the  structure.  Both  of  these  objectives  can  be 
accomplished  by  utilizing  the  results  obtained  from  the  analyses  of  the 
simplified  model. 

In  principle,  the  method  presented  is  similar  to  that  discussed  in 
references  3  and  4.  In  the  present  investigation,  however,  emphasis  is 
placed  on  the  determination  of  the  predominantly  chordwise  modes  and 
frequencies  of  the  structure  in  addition  to  the  modes  and  frequencies 
of  the  overall  beam  motions.  The  results  obtained  from  the  analysis  for 
the  model  with  springs  can  be  used  to  give  both  the  chordwise  and  the 
spanwlse  coupled  frequencies  provided,  of  course,  that  the  appropriate 
choice  of  the  spring  constant  is  made.  The  \incoupled  cross-sectional 
modes  and  frequencies  can  be  obtained  from  the  analysis  by  letting  the 
spring  constant  vanish. 

Since  the  rigid-body  motion  of  the  idealized  structure  depends  only 
on  the  spring  constant  while  the  relative  chordwise  motion  depends  on 
the  flexibility  of  its  members,  the  value  of  the  spring  constant  should 
be  chosen  so  that  the  uncoupled  rigid-body  frequency  (either  transla¬ 
tional  or  torsional)  of  the  idealized  structure  is  the  same  as  the  span- 
wise,,  frequency  of  the  prototype.  For  the  symmetrical  chordwise  modes 
the  spring  constant  should  be  determined  from  the  frequencies  of  the 
spanwlse  bending  modes  by  means  of  the  equation 

k=a>^,BP'tc  (W) 


while  for  the  antlsymoetrical  chordwise  modes  the  spring  constant  should 
be  determined  from  the  spanwlse  torsional  frequencies  by  means  of  the 
equation 


k 


0). 


'r,T 


ptc 

~ 


(^) 
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The  spanwise  frequency  of  the  prototype  should  he  obtained  by  the  most 
appropriate  method  available,  with  the  provision  that  the .method  be  based 
on  the  assumption  that  the  cross  section  of  the  structure  remain 
undistorted. 

At  this  time  a  note  of  caution  as  to  the  application  of  the  method 
presented  is  in  order.  The  use  of  the  idealized  structure  embodies  the 
assumption  that  all  cross  sections  distort  identically.  Hence  the  influ¬ 
ence  of  any  bulkhead,  rib,  or  other  chordwise  restraint  is  not  included. 
Consequently,  the  method  is  applicable  only  to  wings  in  which  the  bulk¬ 
head  spacing  is  not  very  small.  L 

1 

In  order  to  demonstrate  the  use  of  the  method,  both  the  coupled  4 

chordwise  and  coupled  spanwise  modes  and  frequencies  of  a  circular-arc  4 

monocoque  beam  used  in  the  experimental  investigation  of  reference  2  will  4 

be  calculated.  For  the  sake  of  completeness  a  short  description  of  the 
specimen  used  and  the  results  obtained  will  be  presented  in  this  paper. 

The  aluminum  monocoque  beam  was  l40  inches  long  and  4o  inches  wide. 

The  nominal  dimensions  and  average  radius  of  curvature  are  shown  in 

h 

figure  9-  The  beam  had  a  ratio  -2.  of  approximately  0.04  and  a  ratio 

t  2  " 

j  of  0.06735*  The  cross  sections  were  maintained  by  steel  bulk¬ 
heads  located  10  inches  from  each  spanwise  tip.  Eighteen  bending  and 
torsional  modes  were  found,  and  in  addition  a  number  of  so-called  "panel" 
modes  were  located.  (See  figs.  5  and  7  of  ref.  2.) 

From  a  comparison  of  the  spanwise  node  lines  shown  in  reference  2 
and  the  mode  shapes  shown  in  figure  6,  it  can  be  concluded  that  many  of 
these  panel  modes  can  be  described  as  chordwise  modes.  The  variation 
in  position  of  the  chordwise  node  lines  is  the  result  of  coupling  of 
the  symmetrical  chordwise  modes  with  the  beam  bending  modes  or  coupling 
of  the  antisymmetrlcal  chordwise  modes  with  the  beam  torsional  modes. 

Consequently,  a  family  of  chordwise  modes  is  associated  with  each  of 
the  spanwise  modes.  An  additional  node  line  appears  at  each  spanwise 
end  of  some  of  the  chordwise  modes.  These  node  lines  result  from  the 
rigid  tip  bulkheads  used  in  the  construction  of  the  specimens. 

Since  the  covers  of  the  monocoque  beam  were  uniform  in  thickness 
and  the  camber  of  the  cross  section  can  be  closely  approximated  by  a 
parabola,  the  exact  solution  previously  discussed  is  applicable.  Hence 
the  predominantly  chordwise  and  predominantly  spanwise  frequencies  of 
the  monocoque  Including  the  influence  of  coupling  can  be  obtained 
directly  from  figures  5  and  8. 
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Before  this  can  be  done,  however,  the  uncoupled  beam  frequencies 
cUj.  g  and  ^  and  the  uncoupled  cross-sectional  frequency  2. 

be  obtained.  The  uncoupled  chordwise  frequency  co  -i  was  calculated 

X 

with  the  use  of  the  curves  in  figure  ^(a)  or  equation  (39)  foi"  symmetrical 
frequencies  and  equation  (3^)  (with  k  =  O)  for  the  antisymmetrical  case. 
The  uncoupled  chordwise  frequencies  are  listed  at  the  top  of  tables  I 
and  II.  The  bending  frequencies  of  the  beam  were  calculated  in  refer¬ 
ence  2  by  methods  which  include  effects  of  transverse  shear  and  shear 
lag,  and  the  torsional  beam  frequencies  by  methods  which  include  restraint 
of  warping.  These  beam  frequencies  are  listed  in  column  of  tables  I 
and  II. 


In  the  first  column  of  tables  I  and  II  is  shown  a  sketch  of  the 
nodal  pattern  \xrider  consideration.  The  modes  are  grouped  in  pairs, 
except  for  the  first,  with  each  pair  consisting  of  a  spanwise  mode  and 
the  corresponding  first  coupled  chordwise  mode.  (The  higher  chordwise 
modes  were  not  obtained  experimentaliy  and  thus  are  not  included. )  The 


first  mode  in  column  of  each  table  is  the  chordwise  mode  associated 

“cl  “cl 

with  each  of  the  rigid- body  spanwise  modes.  The  ratios  ■  and  ■  ■ 


are  shown  in  column  of  tables  I  and  II, 

spending  frequency  ratios  - -  and  - — — 

“r,B  “r,T 


'r,B 


r,  T 


respectively.  The  corre- 
were  obtained  from  figures  5 


and  8  and  are  tabulated  in  column  (4) .  For  the  predominantly  beam,  modes 
this  ratio  represents  the  reduction  in  beam  frequency  due  to  cross- 
sectional  deformations.  In  an  attempt  to  evaluate  the  results,  the  cal¬ 
culated  coupled  frequency  o)  is  presented  in  column  and  the  experi¬ 
mental  frequency  in  column  (^. 


In  table  I  (for  bending)  the  predominantly  beam  modes  show  little 
effect  due  to  coupling  at  the  low  frequencies.  At  the  higher  modes, 
however,  the  coupling  effect  becomes  appreciable,  and  it  is  encouraging 
to  note  the  good  agreement  between  calculated  and  experimental  frequencies. 
This  good  agreement  at  the  higher  moles  can  also  be  looked  upon  as  an 
indication  that  the  method  is  applicable  to  low- aspect- ratio  wings. 


In  table  II  (for  torsion)  the  predominantly  beam  modes  are  greatly 
affected  by  the  coupling  due  to  cross-sectional  deformations,  even  at 
the  low  frequencies.  However,  the  accuracy  of  prediction  of  the  coupled 
frequencies  at  the  higher  modes  indicates  that  the  applicability  of  the 
method  in  torsion  to  low-aspect- ratio  wings  has  limitations. 


The  calculated  coupled  chordwise  frequencies  show  fair  agreement 
with  experimental  frequencies  and  Indicate  that  a  large  number  of  chord- 
wise  modes  are  present  over  a  narrow  frequency  band,  sometimes  falling 
within  a  few  cycles  of  each  other.  Unfortunately,  as  the  investigation 
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presented  in  reference  2  was  primarily  concerned  with  beam  modes,  the 
experimental  frequencies  of  the  chordwise  modes  were  not  always  obtained 
and  some  are  missing  in  tables  1  and  II.  Further,  because  of  the  com¬ 
plexity  of  the  nodal  patterns,  it  was  difficult  in  some  instances  to 
correlate  the  experimental  and  theoretical  mode  shapes.  Though  the  cor¬ 
relation  between  the  calculated  and  experimental  chordwise  frequencies 
is  good,  the  application  to  low-aspect- ratio  wings  is  open  to  question, 
as  the  accuracy  of  prediction  is  not  well  established  at  the  higher 
frequencies. 

It  is  to  be  noted  that  a  second,  thicker  skinned  monocoque  beam  was  L 

discussed  in  reference  2,  but  for  a  number 

sented  or  discussed  here.  First,  the  f  ^ 

thicker  monocoque  is  felt  to  be  straining  the  limits  of  some  of  the  basic 
assumptions  of  the  present  analysis.  Second,  because  of  a  much  higher 
uncoupled  chordwise  frequency,  the  coupling  effects  on  the  beam  frequencies 
in  bending  would  be  much  smaller  than  in  the  lighter  monocoque.  Further¬ 
more,  the  chord  length  and  radius  of  curvature  used  in  the  calculations 
of  the  uncoupled  beam  frequencies  <1^^, B  ^r  T  questionable, 

as  the  frequencies  of  the  rigid- body  beam  modes  B  calculated  from 

these  dimensions  are  lower  than  those  observed  experimentally.  Since 
the  theories  used  in  these  calculations  do  not  include  all  the  flexibil¬ 
ities  of  the  actual  structure,  the  calculated  values  should  be  higher 
than  the  experimental  values,  rather  than  lower.  It  may  also  be  noted 
thax  the  dimensions  of  maximum  depth,  chord  length,  and  radius  of  curva¬ 
ture  for  the  thicker  beam  are  not  compatible  with  those  for  a  circular 
arc.  When  the  radius  of  curvature  was  altered  by  6  percent  in  order  to 
describe  a  circular  arc,  a  12-percent  change  was  noted  in  the  third  and 
fifth  rigid- body  modes.  More  them  likely,  then,  the  calculated  rigid- 
body  frequencies  of  the  thicker  beam  are  somewhat  higher  in  magnitude 
them  those  calculated  in  reference  2. 


of  reasons  it  has  not  been  pre-  1 

2  ^ 

ratio  of  O.356  for  the  ^ 
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CONCLUDING  REMARKS 


The  vibrational  behavior  of  thin  wings  with  airfoil- type  cross 
sections  has  been  Investigated  by  means  of  an  analysis  of  an  idealized 
structure  consisting  of  two  cunred  beams  mounted  on  springs.  The  results 
obtained  from  an  emalysls  of  this  model  can  be  used  not  only  to  obtain 
the  chordwise  modes  and  frequencies  of  the  prototype  but  also  to  obtain 
an  estimate  of  the  reduction  in  spanwlse  frequency  due  to  the  coupling 
effect  between  cross-sectional  and  spanwlse  motions.  Both  uncoupled 
chordwise  frequencies  and  the  resulting  coupled  chordwise  and  spanwlse 
frequencies  are  presented  for  beams  with  uniform  cover  thickness  and 
parabolic  cross  sections.  A  comparison  is  made  with  existing  experimental 


-P-  -r  4^  H  f 
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data  on  a  circular-arc  monocoque.  The  comparison  shows  good  agreement 
for  the  coupled  beam  frequencies  and  fair  agreement  for  the  coupled 
chordwise  frequencies  at  moderate  aspect  ratios.  For  beams  of  arbitrary- 
cross  section  an  approximate  matrix  procedure  is  also  presented. 


Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Air  Force  Base,  Va. ,  September  21,  I96I. 


22 


APPENDIX 


APPROXIMATE  SOLUTION  OF  AN  ARBITRARY  AIRFOIL  SECTION 


In  this  appendix  an  approximate  solution  for  an  idealized  structure 
of  arbitrary  shape  and  variable  thickness  is  presented.  The  energies  of 
the  system  are  expressed  in  discrete  form  by  writing  all  derivatives  as 
differences  and  all  integrals  as  stimmations  according  to  the  trapezoidal 
rule.  The  deflections  at  discrete  points  are  considered  as  the  unknowns  L 

of  the  system.  After  the  boundary  conditions  are  properly  accounted  for;  ■ 
a  system  of  linear  simultaneous  equations  in  terms  of  these  unknowns  is  I4. 

obtained  by  minimizing  the  energy  function  U  -  T.  4 

4 

From  equations  (4),  (6),  and  (7)  the  expression  for  U  -  T  can  be 
written  as 


The  first  step  in  the  procedure  is  to  choose  equally  spaced  stations 
along  the  x-axis,  numbering  them  from  n  =  0  to  n  =  N.  Integrals  can 
then  be  approximated;  either  in  terms  of  values  of  the  Integrand  at  the 
full  stations 


jT  f(x)dx  =  a(^  +  f^  +  .  .  .  +  fjj_^  +  (A2) 

or  in  terms  of  the  values  at  the  half  stations 


Jq  -  ^[^1/2  ^5/2  +  •  •  •  +  ^N-(l/2)] 

where  A  is  the  equal  distsince  between  stations. 

The  first  Integral  in  equation  (Al)  is  approximated  by  means  of 
equation  (A2),  as  the  values  of  the  second  derivative  appearing  in  the 
integrand  are  best  approximated  by  differences  at  the  full  stations: 


H  -d"  -d" 


1 
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w 


w  ,  T  -  2w  +  V  T 

n+1  n  n-1 


n 


(A4) 


Consequently,  the  first  integral  can  be  vritten  in  matrix  form  as  follows: 


El(v")^dx  =  ^[w J  ^2*]  '  [eiJ  [02' 


w 


(A5) 


where  the  definitions  of  the  matrices  are  given  below,  and  the  number  of 
rows  and  columns,  in  that  order,  are  indicated  beside  each: 


l_w  =  j^w_p  Wq,  .  .  .,  (l)x(N+3)  (A6) 


w 


(N+3)x(l)  (AT) 


(N+l)x(N+5)  (a8) 


1-2  1 
1-2  1 


(n+i)x(n+i)  (A9) 


A  prime  on  a  matrix  denotes  the  transpose  of  the  matrix. 


2k 


Note  that  the  values  of  w  at  stations  corresponding  to  n  =  -1 
and  n  =  N  +  1  (points  not  on  the  structure)  appear  in  the  matrices. 
These  points  can  be  evaluated  with  the  use  of  the  finite-difference  form 
of  the  boundary  conditions  in  equations  (lO)  and  (12): 


EIw" 


n=N 


w. 


n+1 


'^N-l 


=  ( 


and 


ETw" 


n=0 


Wj  -  2wo  +  w,^ 

a2 


=  0 


(AlO) 


(All) 


With  the  conditions  expressed  in  equations  (AlO)  and  (All),  equa¬ 
tion  (A5)  can  be  rewritten  as 


L 

1 

4 

k 

4 


(A12) 


where  now 


[wj  =  I^Wq,  Wp 


•  • 


(1)X(N+1) 


(A15) 


1-2  1 


1-2  1 


1-2  1 
1-2  1 


(N-1)x(N+1)  (A14) 


(N-l)x(N-l)  (A15) 
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L 

1 

k 


The  second  integral  in  equation  (Al)  is  approximated  hy  the  use  of 
the  values  at  the  half  stations  (eq.  (Aj)),  as  the  value  of  the  first 
derivative  appearing  in  the  integrand  is,  in  turn,  best  approximated  by 
a  half-station  expression: 


w' 


^  ^n+1  -  ^n 
n+(l/2) ~  A 


(Al6) 


Thus,  by  using  both  the  half-station  expressions  (eqs.  (AJ)  and  (A16)) 
the  second  integral  can  be  written  in  terms  of  the  displacements  at  the 
full  stations: 


i  Eth’w'docJ^  h’w'dx  =  iLwj[Dj'[h’J[Et]^h^[^D^|  w|  (Al?) 


Here 


H  ■ 

Wjj 

j  (1)X(N+1) 

(A18) 

(N+1)X(1) 

(A19) 

-1  1 

-1  1 

-1  1 

-1  1 

-1_ 

(n)x(n+i) 

! 

(A20) 

\/2 


b/2 


^-(5/2) 


%-(l/2) 


(N)x(N)  (A21) 
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<"*>1/2 

— 

OJ 

1 — I 

<"*>3/2 

(Et)^/2  •  •  • 

^®’'^^N-(l/2) 

<"*>H-(l/2)  •  •  ■ 

(N)X(N) 


(A22) 


The  third  term  in  equation  (Al),  or  the  potential  energy  of  the 
spring  mounts,  is  given  in  matrix  form  hy 


i  k^v(c)^  +  i  kgwCo)^ 


w 


where 


(n+i)x(n+i) 


(A25) 


(A24) 


The  last  integral  in  equation  (Al)  is  approximated  hy  equation  (A2): 


2 


where 


(A25) 


(N+1)x(N+1)  (A26) 


i-P  iH  -4"  -4- 


27 


L 

1 

4 
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With  equations  (A12),  (A17),  (A23),  and  (A25)  the  energy  function 
(eq.  (Al))  can  be  written  as 


fr  1  /,^2| 

r 

H 

[[-] -  tI 

L®! 

w 

where 


[b]  =  a[pi 


N  n 
h 

■sJ  L 


Dn 


(A27) 

(A28) 

(A29) 


Application  of  the  variation  condition  expressed  in  equation  (8), 

8(U  -  T)  =  0^  where  the  variations  are  taken  with  respect  to  each  w, 
results  in 

([ 

Equation  (AJO)  can 


:^]  ^  H'l  - 


B 


w 


=  0 


(A30) 


be  rewritten  in  a  more  usable  form  as 


since 


(A51) 


(A32) 


6ind  the  transpose  of  A  is 

In  order  to  establish  the  number  of  stations  necessary  to  achieve 
reasonable  accuracy  with  the  finite-difference  method,  a  brief  comparison 
was  made  between  the  deflection  shape  of  the  first  symmetrical  mode  as 
calculated  from  the  exact  solution  presented  for  the  monocoque  wing  of 
symmetrical  parabolic  cross  section  and  the  shape  as  calculated  by 
finite -difference  methods.  Since  a  symmetrical  mode  was  chosen  for  the 
comparison,  only  half  of  the  doubly  symmetrical  parabolic-arc  cross 
section  was  considered  and  the  origin  was  moved  to  the  midchord.  (The 
limits  of  the  Integrals  in  eq.  (Al)  thus  become  0  and  c/2  instead  of 
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0  and  c. )  With  the  additional  simplification  o?'  constant  cover  thick¬ 
ness  and  constant  material  properties,  equation  (A31)  reduces  to 


with 


|w|  =  — ^[Bj-^2[A]|w| 
cT)pta> 


(n+i)x(n+i) 


(A34) 


(A55) 


4ei 

c5t,5 


+ 


32ho2Et 

c5 


(N+1)x(N+1)  (A56) 


A  =  I  Ti  (A5T) 

Since  the  origin  is  now  located  at  the  midchord  the  boundary  condition 
originally  given  by  equation  (All)  (at  x  =  O)  is  replaced  by  that  of 
symmetry  at  the  midchord: 


L 

1 

k 

k 

k 


w_i  = 


(A38) 


and  D2  in  equation  (A55)  becomes  Dj  in  equation  (A56): 


-2 

1 


2 

-2  1 
1-2  1 


1-2  1 
1-2  1 
1-2  1 


(N)x(N+l)  (A59) 
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Furthermore, 


1 

1 


(N)X(N) 


2 


3 


(N)X(N) 


1 


1  1 


(N)X(N) 


(A40) 


(A4l). 


(A42) 


-1  1 
-1  1 
-1  1 


(n)x(n+i)  (a43) 


(N+l)x(N+l)  (a44) 
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In  figure  10  a  plot  of  the  first  symmetrical  mode  shape  of  the 

2 

parabolic-arc  monocoque  for  k  =  0  and  f I  =  0.0i8  is  shown.  The 

\^o/ 

curve  was  calculated  from  equations  (59)  (4o)  of  the  exact  solution, 

and  points  were  obtained  from  equation  (a45)  of  the  numerical  procedure 
for  6,  S,  and  11  stations  across  the  half-chord.  The  error  involved  in 
predicting  the  overall  frequencies,  in  comparison  with  the  exact  solu¬ 
tion,  is  6  percent  for  6  stations,  2.5  percent  for  9  stations,  and  less 
than  2  percent  for  11  stations  across  the  half-chord. 
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Figure  10.-  Comparison  of  finite-difference  method  with  exact  method  for  the  first  symmetrical 

mode  shape,  k  =  Oj  (•r~\  =  0.018. 
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